A model for elasticity and pinning of domain walls in helical magnets 
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The theory of elasticity and coercivity of domain walls in both, centrosymmetric and non- 
centrosymmetric helical magnets, is developed. Generically these walls consist of regular arrays 
of magnetic vortex lines. Exceptions are walls oriented along high symmetry directions. The elas- 
ticity of the latter turns out to be non-local below their roughening transition, pinning of these walls 
by disorder is negligible. In contrast, weak anisotropy, breaking the f/(l) symmetry of the helix, 
results in domain wall pinning by the bulk. Walls of other orientations consist of pairs of vortex 
lines separated by vortex free parts. These are strongly pinned by disorder. The application to 
chiral liquid crystals is briefly discussed. 
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Introduction. — Pinning plays a key role in condensed 
matter systems: it restores the state of zero resistance in 
type-II superconductors by anchoring flux lines [T] and 
hardens steel by blocking the motion of dislocations, but 
it also prevents charge density waves to become ideal 
conductors [2]- In ferroelectrics and ferromagnets pin- 
ning of domain walls (DWs) influences their hardness and 
switching behavior [31 11] , strongly relevant for potential 
applications as storage media [S]. In all cases pinning 
and consequently coercivity result from the competition 
of the impurity potential, which favors deformations of 
the condensed structure, and the rigidity of the latter, 
which penalize them. The appearance of a non-zero co- 
ercive force requires the emergence of multistability of 
the resulting effective potential |B]. Its detailed descrip- 
tion demands therefore the use of elaborate techniques 
like the functional renormalization group (FRG) method 
[7] . Recently a new type of DWs in helical magnets has 
been predicted |5]. Helical magnets are abundant, occur- 
ring as metals and alloys [5HT3]. semiconductors [13] 
and multiferroics [T2H22]- The latter group is most in- 
teresting for applications [16]. It was shown in |8] that 
for almost all orientations DWs in these systems consist 
of a regular array of magnetic vortex lines. These walls 
can be driven by currents and in multiferroics by elec- 
tric fields [S]- Vortex walls have indeed be seen in circu- 
larly polarized X-rays in Ho [10] , and by Lorentz-TEM in 
FeGe [13]. Only for special orientations DWs are vortex 
free. The latter have been studied previously by Hubert 
[3]. However, at least in films, Hubert walls (HWs) were 
not found [23]. Finally we mention that there are close 
analogies between the ordering of helimagnets and of chi- 
ral nematic and smectic liquid crystals [24l [25] such that 
DWs with similar properties have to be expected there 
as wells well . 

Outline and results. — In the present letter we inves- 
tigate the elasticity and impurity pinning of DWs in 
helimagnets. We show that the standard model for 
(centrosymmetric) helimagnets can be considerably 
simplified by writing it in a form resembling the London 
theory of superconductors. In this description magnetic 



vortices result from a fictitious magnetic field acting only 
inside the DW. HWs are found to exhibit non-local elas- 
ticity below the roughening transition temperature Tr. 
They are therefore too stiff to be pinned by impurities. 
Thus metastable HWs can disappear easily from the 
sample, which may explain why HWs were not observed 
so far in experiment. Roughening of HWs occurs by 
entropy driven proliferation of steps consisting of pairs 
of vortex lines of the same vorticity. Against expectation 
these vortices neither repel each other nor do the have an 
energy increasing logarithmically with system size. An 
additional disorder driven roughening of HWs is shown 
to have a weak effect only. On the contrary, vortex 
walls exhibit local elasticity and are strongly pinned 
by the local coupling to the disorder. This is shown by 
using a FRG method. The estimated size of metastable 
domains is in fair agreement with the experiment. We 
finally discuss a non-local (bulk) pinning mechanism for 
HWs which exists in the presence of a weak anisotropy, 
breaking the U{1) symmetry of the magnetic structure. 
Although most of the derivations are presented in detail 
for centrosymmetric systems, the results transfer to the 
non-centrosymmetric case as well. Finally we discuss 
some conclusions for DWs in chiral liquid crystals. 
Hamiltonian. — To describe helical magnets we use the 
appropriate Ginzburg-Landau-Hamiltonian 'H[m(r)]. 
m — {mx,my,mz) denotes the magnetization. Below 
we will parametrize m by the angle ^, rux + imy — e*"^, 
assuming = 1 and = 0. Since in helimagnets 
both time and space inversion symmetry are broken, 
their paraphase can be centrosymmetric or non- 
centrosymmetric. Centrosymmetry requires invariance 
with respect to space and time inversion, i.e. r — > — r 
and m — >■ — m. If there are only two modulation vectors 
q = ±{9/a)x, as in most centrosymmetric systems, one 
finds up to quadratic terms in m [3 
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Here /j. = / d'^r and Vj^ — y dy + z dz. 9 and a denote 
the angle between adjacent spins along the x-direction 
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and the lattice constant. The continuum approach 
is vahd provided 9 <^ n. Experimentally one finds 
0.27 — 0.73 under ambient conditions [9] and 9-^0 
under uniaxial pressure |26) . In systems where indirect 
RKKY exchange between 4f electrons [23 results in 
nearest neighbor ferromagnetic (J > 0) and next 
nearest neighbor anti- ferromagnetic ( J' < 0) interaction, 
9 = arccos( J/4| J'l) in the ground state. 
With the replacement nix + "ifny — e"^, the Hamiltonian 
^ includes a non-linear term [{dx(f>)'^ — q^]'^, making 
calculations cumbersome. We will therefore resort in 
the following to an approximation and replace this 
expression by Aq^ldxcj) — A(r)]^. ^(r) is assumed to 
be constant inside a domain {A — ±q) and to change 
smoothly from —q to g on a scale q~^ when crossing the 
DW parallel to its normal h. The resulting Hamiltonian 
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resembles the London theory of type-II superconductors. 
A = (A(r),0, 0) plays the role of a vector potential 
which generates a Active magnetic field B = V x A 
acting only inside the DW where it creates vortices, q^^ 
corresponds to the London penetration lengths. 
Elasticity of HWs. — Calculation of pinning forces re- 
quires knowledge of the elasticity of the DW [B]. The 
energy of long wave length elastic DW distortions u from 
a planar reference configuration can be written as 



OuiCHO- (3) 



^ = (Ci, C2) are the coordinates perpendicular to the DW 
normal n. To determine GiC) we consider first HWs, 
assuming that their distortions are free of vortices (but 
see below). In this case the saddle point equation for 0, 
following from ([T]), can be solved exactly [3j. For an 
isolated planar wall the solution is 00 (a;) = \ucosh{qx). 
The choice ^(r) = qtanh{qnr) in ([2| reproduces 0o(x). 
The solution of the saddle point equation, corresponding 
to ([2| , for a distorted HW can be found using the Ansatz 
invented in |28j for contact lines in wetting problems 



dx'A{x') + 
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Here = J (P k / {2tt)'^ . The first term describes the 
unperturbed phase field, the second term its corrections 
up to a distance |k|^^ from the average wall position, ak 
is to be determined from (j){u{^), ^) = 0. To lowest order 
in u, Qfk = — Ukg, provided |k| < q. Plugging Q back 
into ([2| we get for the Fourier transform of GiC) 



g^i(k)« J|g|3(4|k|+2akV3) 



(5) 



([5]) is valid for distortions \qu\ > 1. The dominant term 
~ |k| results from the long range interaction of the wall 
which a makes the DW very stiff and suppresses pinning 



of HWs, as we will show below. The second term in ([s]) 
is the contribution from the increased surface area due to 
the wall distortion and is relevant only for large k. In real 
space the energy expression ([s]) for the HW is non-local, 
Q^^[r) ~ J(|q|/|r|)'^. This is in strong contrast to Bloch 
and Neel walls, whose elasticity is related to anisotropy 
and strictly local. 

Roughening transition of HWs. — Next we consider the 
possibility of a roughening transition of the HW which 
would render elasticity short range |29j . To this aim we 
consider a HW with a step at ^1 = 0, parallel to the (2- 
axis [see Figjl] right]. It can be described by the function 

0s (a;, Ci ) = 00 {x-TTT/q)+ TTT sigu {x - Tir/q) , (6) 

where r — signal. (j)s{x,C,i) is smooth across Ci = for 
|a;| 3> Tr/g. On the contrary, in the region |a;| < rr/q the 
sign of dx<ps is opposite for Ci ^ 0. Hence the integral 
along a contour C inclosing the step gives §Q<t> = 47r, i.e. 
the step consists of two vortices. Since in this construc- 
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FIG. 1. Cross section of a HWs (dashed line). Left panel: 
planar wall, right panel: wall with step before relaxation. 
The arrows denote the orientation of m. For systems where 
m is confined to the plane perpendicular to x, m has been 
rotated by 7r/2 around for better visibility. The contour C 
encloses two vortices. 

tion the vortex configuration is restricted to a narrow 
slice of width a, the step energy per unit length, eji, is of 
the order J/{a9). Further relaxation of the configuration 
0s by allowing the vortex to extend over a larger region 
can only decrease Eh- A variational calculation gives then 
£h ~ (7r%/5/4) Jk/o where = ln(7r/26'). The fact that 
Eh is finite, as well as the attraction of the two vortices 
are surprising: typically vortices repel each other and are 
associated with an energy increasing logarithmically with 
the system size. In the present case however this energy 
is saved by the presence of the HW. For energetic reasons 
the two vortices also do not separate. 
Using the step energy eh in the results for the roughening 
transition in the ASOS model [301 one obtains for the 
roughening transition temperature of the HW 



2.2jK/kB ~ kTjs! . 



(7) 



In the second step T/j is compared with the Neel temper- 
ature Tm- For T/v we adopt here the result for the 3D 
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XY-model, T/v ~ 2.2J/kB [31], corresponding to the use 
of ([2]). This approximation is restricted to a region not 
too close to the Lifshitz point 9 = 0. Above Tr the HW 
exhibits short range elasticity. However, since Tr is of 
the order T/v, the non-local elasticity dominates over a 
large temperature region. 

Elasticity of vortex walls. — Next we consider a DW of 
arbitrary orientation, which include vortices. Denoting 
by a = arccos(ri- x) the angle between the x-axis and the 
wall normal h, for not too large a, the wall consists of a 
finite density of those steps considered before. Its surface 
tension a{a) can then be written as 

<7{ce) ~ [an + £int(a/ tana)] | cosaj + ay\ sinaj. (8) 

an = 2J9^/{3a'^) and = e//0/(27ra) are the vortex 
tension of the Hubert and the pure vortex wall {fix = 0), 
respectively, eint(^) describes the step interaction. At 
T = eint(^) ~ cTv exp {—q£), whereas at T > eint(^) ^ 
T/P due to collisions of meandering steps [^ . 
To find the elastic properties of the tilted wall we choose 
Ci = C2 X C2 II X x h, such that the vortex lines are 
parallel to the ^2 direction. For an infinitesimal homoge- 
neous distortion di^^u = e one obtains for the change of 
the surface energy density 

- a(a) « a'{a)€ + ^ [a{a) + o{a)"] . (9) 

The linear term in e vanishes at finite temperatures since 
Tr = for the vicinal surface considered here. This 
follows from the fact that HWs are structureless in the 
plane orthogonal to x and hence steps can meander freely, 
leading to a rough surface [53] . A distortion d^^ u leads to 
an expression without derivative terms since the surface 
tension depends only on a. The total elastic energy for 
the vortex wall can therefore be written as 
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where JCi(a)/a^ = a{a)+a" {a) and JC2(a)/a^ = a{a). 
Thus, vortex walls exhibit conventional elasticity. The 
mixed term in the elastic energy vanishes for symmetry 
reasons. Note that Ci only depends on the vortex in- 
teraction which is small for small a. A corresponding 
calculation for a w 7r/2 is more difficult. The height 
h of steps in pure vortex walls can take any value, but 
steps cannot meander freely since vortex walls have a 
structure periodic in the ^-direction. However, since the 
step energy e^(/i) is small, vortices can almost freely slide 
against each other; hence Tr is low in this case [32]. 
Pinning of HWs. — DWs can be pinned by impurities. 
The statistical pinning theory of DWs with local elastic- 
ity has been developed some time ago [7] . In the present 
context it applies to vortex walls. We will therefore con- 
centrate in the following on the new features arising for 
pinning of HWs below their roughening transition. We 



assume the presence of non-magnetic impurities of con- 
centration riimp which dilute the system and hence con- 
tribute a term 'Himp — —Jl X^i ^[4>{''^i^ "")] to the energy. 
J 8 /a^ is the energy density of the domain wall, 7 and 
denote the impurity strength and position, respectively. 
For HWs £h « a^Wo{x - u)f 12 [3]. The local pinning 
force density follows then from 
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J2SiC-(:^)Hx^-u) (11) 



where = {xi,Ci) and J-{x) = ,770* [(/)q^(x)] '/2. Aver- 
aging over the impurity positions we get (/) = and 

(/(C, n)f{C, u')) = S{C - C')Ao(« - u'). (12) 

Here (...) denotes the disorder average and Ao(m) = 
TT-imp J dxT{x — u)J-(x). Second order perturbation the- 
ory gives for the pinning threshold of a driven DW [7] 



/pin --A^(U^ ±0)^(0). 



(13) 



Here ± sign denotes the sign of the driving force. Gener- 

ically A'oiO) / dx{T^{x)Y = for analytic A(u), 

since T{x — t- ±00) = 0. Then /pin = 0, i.e. there is 
no coercivity from perturbation theory. Instead we can 
use a FRG approach for A(u) , provided the dimension 
D{= 2) of the DW is below the critical dimension Dc- 
The resulting pinning strength then increases exponen- 
tially with Dc — D. Dc follows from the comparison of 
the energy Eg\ of an elastic strain l/{qL), considered on 
a sufficiently large scale L, with the fluctuation of the 
pinning energy Epi^ ^ L^^^ on this scale. If Egi < i?pi„, 
the DW can adapt to the disorder and hence accomm- 
modate to a potential valley where it gets pinned. In the 
opposite case the DW is to stiff to stay in one valley. By 
crossing the rugged energy landscape, potential forces on 
the DW show either sign such that the resulting pinning 
force is proportional to the square root of its area. It is 
therefore surpassed by the driving force, which is typi- 
cally proportional to the area. For systems with short 
range elasticity Ec\ ^ L^^^, and hence Ec\ < -Epin for 
D < Dc = 4. On the contrary, for HWs E^i ^ L^^^ and 
hence Dc = 2, resulting in exponentially weak pinning, 
as we will show now. 

Following the calculation scheme used in [^ one can show 
that the effective force correlator Ai(u) on scale L = ae^ 
obeys the RG flow equation (c = l/[64:nJ^q^]) 

= c^.Mu) [2A,(0) - Mu)] . (14) 

Introducing the curvature of the force correlator at u = 
0, gi = A^'(O), we find from ^ and A'(0) = that 
g obeys the equation dgi/di 



Qcgf. Integration with 
the generic initial value 50 < shows, that g develops a 
pole on scale C = aexp[l/(6c5o)]- On larger scale Af(M) 
exhibits a cusp. In this region ([l4| can be solved with the 
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ansatz Ae{u) = c-^A'^£-'^+'^''A*{uA-^e-t') which gives 

^uA*' - A* " 



for A*(u) the relation 
(1-2/i- A*")A 



A* 



0. (15) 



We have chosen A such that A*(0) = 1. /i = 1/3 can be 
found from the fact that de J du/S.i{u) = 0. For small u 
one then gets A*(u) = - \u\/y/3 + 2u'^ /9. Thus A*(±0) - 
— signw and oo > g* > 0. We can now apply (131 using 
the renormalized function A.£(u) on scales larger C. The 
latter plays the role of the short length scale cut-off of 
the renormalized theory. With go «— 2.4ci, 
one obtains finally for the coercive force 



-'imp ^ 



\j pi 



(H) 
pin 



^^^^exp[-ci/(|^|^Ci„,p)]. (16) 



Cimp — 7^"-impa'^ is proportional to the volume fraction 
of the impurities and ci ~ 14.2. Since both 9, Cimp <C 1, 
pinning of HWs by direct interaction with impurities is 
completely negligible, as long as one is below Tfj. This is 
a direct consequence of their non-local elasticity. It im- 
plies that after a quench to a metastable multi-domain 
state HWs will quickly disappear from the sample. In- 
deed, in films of helical magnets with film plane perpen- 
dicular to the helical axis, domains were found to extend 
over the whole film width [53] . 

We have also studied the possibility of disorder driven 
roughening of HWs which would render its elasticity lo- 
cal [33]. Using arguments similar to those used in [34] 
we found that the step energy surrounding an terrace of 
linear size L, is now lowered due the disorder induced cor- 



rugation by Eh, imp - £h ^ -{J/ a 



--imp 



ln(L/a))2/3. 



SH.imp vanishes on scales L > aexp[£/fa/( J|6'|^Cimp;j. 
Hence terraces are generated spontaneously by the disor- 
der. Because of the exponentially large length scale this 
effect will hardly be seen and hence HWs remain flat and 
elasticity non-local as long as T < Tr. 
Bulk pinning of HWs. — In systems where the U{1) sym- 
metry of the magnetic structure is broken by weak 
anisotropy there is also a non-local interaction between 
impurities and the HW, as we will show now. Moving a 
rigid HW requires the rotation of all spins in at least one 
half space [see Figj2] . Such a rotation does not cost en- 
ergy even in the presence of impurities, as long as the sys- 
tem has U{1) symmetry and the interaction with impuri- 
ties does not depend on the spin direction, as we assume 
here. Most of the experimental systems have however a 
weak anisotropy of the form 

ny = - j{v/a^) cosipcj)), w > 0, (17) 

which generates a weak modulation of the wave vector 

V 

q{x) « g [1 - iysmp{qx + , - ^ < 1. (18) 

This results in turn in a modulation of the energy density 
of the helix Sh ^ Ja{q'{x))^ . Dilution by non-magnetic 
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FIG. 2. Left: Spatial variation of the top of magnetization 
vector (black) along the a;-axis with one HW (kink), separat- 
ing a left-handed from a right-handed helix. Right: Schematic 
plot of the phase profile (ftix) in the case of (strong) bulk pin- 
ning for non-zero anisotropy v. The Xi denote the impurity 
positions. The dashed line shows (f>{x) for v = Q. 



impurities acts now strongest close to the maximum of 
|g'(a;)|, i.e. at qx = 4>q-\- ni: /p. Depending on the sign of 
the interaction, these regions are either attracted or re- 
pelled by impurities. Since impurities are randomly dis- 
tributed, the energy of a rigid helix will again not depend 
on a;o and can hence be rotated without changing its en- 
ergy [see dashed line in Figj2] right panel]. The situation 
changes if the helix can deform elastically by changing 
its wave vector by 5q{v). A stretch of the helix by 5q on 
scale L costs an energy ~ JL^O^Sq'^ /a. But since impu- 
rities are distributed randomly, the typical energy gain 
due to the stretch is of the order Jv(\l^^{L/a)'^/'^ . A dis- 
tortion of the order 5q « tt/L is therefore favorable on 
the scales larger than « a(7r^^)^/(^;^Cimp). On these 
scales the phase profile 0(x) has adapted to the disor- 
der and resists now a change of the phase 0o- This can 
most easily be seen in the case of an infinitely soft one- 
dimensional system, corresponding to a strong pinning 
situation. There the condition of maximal pinning energy 
gain, 0(xi)-|-0o = 27rni/p, can be fulfilled simultaneously 
at all impurity positions Xi [see solid line in Figj2j right 
panel] . The integers are chosen to minimize the resid- 
ual elastic energy. A global spin rotation, corresponding 
to a shift of the modulated structure, would now change 
the energy. This leads to pinning. The pinning force 
per unit volume, resisting rotation of the helix, is then 

and the coercive force 



(bulk) 



(H.huWi) 
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,^./(^V|0|^). 



(19) 



Here is the typical distance between HWs. 
Pinning of vortex walls. — Their elasticity is short range, 
shifting Dc to 4. The relation between the Larkin length 
and AJj'(O) follows in complete analogy to the calculation 
for the HW if we replace ([s]) by ( 10 1 in the calculation of 
the flow of Af. This gives Cl « y/CiC2 / {a^\K{^)\) ■ 
To obtain A"(0) we have to determine the energy density 
~ £y of the vortex wall. Since in this case there is no 
exact solution of the saddle point equation, we resort 
again to the variational calculation used in Clearly 
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£v{x,y) of a wall in the x-z plane is periodic in x with 
period tt/q. We therefore can restrict ourselves to the 
region \x\ < TT/{2\q\) where we get (c2 = 5/64) 



£v{x,y) 2c2 [aKx/ {c2y^ 



(20) 



(201 reproduces — ^ £{x,y)q/aTr, as it should be. 



inning force J-{x,hr) — Jj£y{x — nh,nr) and disorder 
correlator A"(x, 0) ^ J^Cimpa/ {xkY are now modulated, 
reflecting the vortex array inside the DW. The main con- 
tribution to the pinning force comes from the center of 
the vortices. Averaging A"(a;,0) over the one vortex pe- 
riod one obtains for the coercive force 



(v) 



K ) (21) 



where C3 ~ 1.7 x 10^'"'. The result can be used to calcu- 
late the typical size (perpendicular to x) of metastable 
domains. Assuming equilibrium between surface tension 
and the pinning force, one obtains from = fp^i^R for 
the size of metastable domains R ~ C4a\9\K^ /cimp. Here 
C4 « 4.6 X 10^. Using ni^p « lO^^m-^, a = 4 x 10-^°m, 
J — 4 and 6 = 7r/6, one obtains R « 200/j,m, in reason- 
able agreement with experiments [TU] . It should be noted 
that in contrast to simpler scaling arguments, the small 
factor C3 can only be found from a FRG. 
N on- centra symmetric case. — In the generic case [S] 

?^ = ^^[(Vm)2 + |q|m(Vxm)]. (22) 

The direction of q = 6q/a is fixed by a weak (cubic) 
anisotropy of the order 9'^. 9{<^ 1) is here proportional to 
the spin-orbit coupling. The anisotropy term can there- 
fore be neglected otherwise. HWs are characterized by a 
normal fi obeying r;,q_|_ = nq- . q+ , q_ are the wave vec- 
tors of the adjacent domains. To be specific we consider 
Qi = (9^7 0, ±gz(x)) with qz{x) changing smoothly from 
qz to —qz over a region of size q~^ when crossing the wall 
[8]. Expressing m — y cos (j) + q{x) x ?/sin0, and ignoring 
terms which are non-zero only inside the HW (and hence 
do not contribute to non-local elasticity), we can rewrite 
( 22 1 in the same form as ^ 



^^[V0-q(x)] 



(23) 



Thus, HWs in non-centrosymmetric systems show long 
range elasticity as well. All other conclusions made for 
HW and vortex walls transfer correspondingly. 
Liquid crystals. — Chiral nematic and smectic phases of 
liquid crystals exhibit helical phase [24l [25]. Their 
description is similar to that used in the non- 



centrosymmetric case, (22), provided m is replaced by 



the director n, and J by the Frank constant. Further, in 
contrast to helimagnets, the direction of q is not fixed in 
space since chirality is introduced through the chirality 



of the molecules. DWs of the type described above will 
occur however as grain boundaries between phases with 
different q direction. A detailed discussion is however 
beyond the scope of the present paper. 
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